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ABSTRACT 

Compact relativistic stars allow us to study the nature of matter under extreme 
conditions, probing regions of parameter space that are otherwise inaccessible. Nu- 
clear theory in this regime is not well constrained: one key issue is whether neutron 
stars are in fact composed primarily of strange quark matter. Distinguishing the two 
possibilities, however, has been difhcult. The recent detection of seismic vibrations in 
the aftermath of giant flares from two magnetars (highly magnetized compact stars) 
is a major breakthrough. The oscillations excited seem likely to involve the stellar 
crust, the properties of which differ dramatically for strange stars. We show that the 
resulting mode frequencies cannot be reconciled with the observations for reasonable 
magnetar parameters. Ruling out strange star models would place a strong constraint 
on models of dense quark matter. 
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The properties of compact stars offer the best oppor- 
tunity to study the phase diagram of Quantum Chromody- 
namics (QCD) at extreme densities, a region that is poorly 
constrained. One hypothesis, based on the possibility that 
strange quark matter (SQM) could be absolutely stable, is 
, that c ompact stars a re made almost entirely of deconfined 

■ SQM (|Witten|[l983 '). In this scenario, terrestrial matter 
] made of nuclei is only metastable, albeit with an essentially 

■ infinite lifetime because it would require mu ltiple weak in - 
' teractions to convert normal matter to SQM (| Weber Il2006h . 

, However, the extreme ambient conditions characteristic of 
compact stars can facilitate the conversion of ordinary mat- 
ter to SQM on short enough time-scales so that some or all 
compact stars could be strange stars. Distinguishing strange 
stars from neutron stars, however, has proved difficult. Sig- 
natures such as potentially smaller radii, and a supposed in- 
ability to glitch, are fraught with serious observational and 
theoretical uncertainty. 

New observations, however, could change this. Cer- 
tain magnetars (compact stars with magnetic fields ^ 10^'' 
G), the Soft Gamma Repeaters (SGRs), exhibit gamma- 
ray fiares powered by field decay ^Duncan & Thompson 

1 19921 ). Timing analysis of rare giant flares, which have peak 
luminosities 10''* — 10*^ erg s~' and decaying tails last- 
ing several minutes, has recently revealed high frequency 
Quasi- Periodic Oscillations (QPOs). The 2004 giant flare 
of SGR 1806-20 shows QPOs at 1 8 26, 30, 92, 150, 625 
Hz and higher |lsrael et al~l l2005l : IWatts fc Strohmaver I 

I2OO6I : IStrohmaver fc Watts 1 120061 ). whilst the 1998 giant 



flare fro m SGR 1900-H4 has QPO s at 28, 53, 84 and 
155 Hz (|Strohmaver fc Watts 1 12005| ). The most promis- 
ing model involves seismic vibrations, triggered by a star- 
quake associated with the giant flare. The lowest fre- 
quency 18, 26 Hz QPOs fit predictions for Alfven modes 
of the core. For the higher frequencies, attention has fo- 
cused on toroidal shear modes of the crust (or their 
global magneto-elastic equivalent). These have frequencies 
that are, for neutron stars, a good match to t he obser- 
vations jM cDermo tt. van Horn fc Hansen I 19881: Duncan 
1 19981 : iPiro 1 i2005i: IGlampedakis. Samuelsson fc Andersson 



200^; ISamuelsson fc Andersson I |2007^ ■ The 28-155 Hz 
QPOs would be n = modes (no radial node) with differing 
angular quantum number I. The 625 Hz QPO is consistent 
with being the n = 1 first radial overtone. 

The likely dependence on crust properties is very 
exciting. Originally, strange stars were expected to be 
devoid of a solid crust and were characterized by an 
ultra-dense quark liquid exten ding up to the surface 
(|Haensel. Zdunik fc Pichon Ill986l ). Such bare strange quark 
stars could not account for torsional shear oscillations since 
there is no solid region in the vicinity of the surface. How- 
ever, strange stars can have solid crusts. One possibility 
l|Alcock. Farhi fc 01into1ll986l ) is that the strange star has 
a thin crust of normal nuclear material extending down to 
neutron drip at density p « 4 x 10^' g cm"'' (for a neutron 
star the crust extends beyond neutron drip to p « lO'* g 
cm~^), suspended above liqui d SQM by an enormous elec- 
tric field. A more recent model (jjaikumar. Reddv fc Steiner I 
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I2OO6I ) posits a crust in which nuggets of SQM are embed- 
ded in a uniform electron background. The strange star 
crusts have different shear speeds and are thinner than neu- 
tron star crusts: both factors affect shear mode frequencies, 
crust thickness AR being of particular importance to the 
radial overtones (lianscn & Cioffi 198(1). Magnetar seismol- 
ogy may therefore offer a robust means of distinguishing 
strange stars from neutron stars. 



1 MODEL 

Most previous studies computing torsional shear modes of 
the stellar crust assume free slip over the fluid core. But 
if the strong magnetic field couples crust and core together 
then one should instead consider global magneto-elastic per- 
turbations. Recent work on the coupled problem, however, 
suggests that the vibrations most likely to be observed 
have very similar frequencies to those computed in the un- 
coup l ed problem jGlampedakis. Samuelsson fc Andersson I 
|2006I :Il evm 1l2007l n The natural crust frequencies are there- 
fore physically relevant, so for simplicity we n eglect coupling 
and apply zero traction boundary conditions ijCarroll et al. I 

1 19861 ). 

FoUowing lPirol Hooi) we use a plane-parallel geometry 
with constant gravitational acceleration g = GM/R^ in the 
vertical direction z, M and R being the mass and radius 
star. The Newtonian equations of hydrostatic equilibrium 
determine the crust density profile. Using a slab rather than 
spherical geometry allows us to incorporate the magnetic 
field without difhculty; we assume a constant field B = Bz. 
For pure toroidal shear modes, which are incompressible and 
have no vertical component of displacement, the perturba- 
tion equations for the horizontal displacement ^ reduce to 



) pR2 



? = (1) 



where fi is the shear modulus and va ~ -B/(47rp)^''^ 
the Alfven speed. Shear speed Vs = (fi/p)^^'^. We have 
assumed a periodic time dependence exp(jaji), lj being 
the frequency. We subsequently correct for gravitational 
redshift to obtain the observed frequency. Primes indi- 
cate derivatives with respect to 2. We have used VxC — 
— [{I + 2){l — 1)/R'^]^ to mimic a spherical geometry, I 
being the standard angular quantum numbe r. The scal- 
ing with I differs from that used by IPiro I (|2005l ) but 
gives better agreement between slab and spherical mod- 
els in the zer o field limi t ('McDermot t. van Horn fc Hansen I 
Li988l Samu elsson fc A n dersson ,,2007^ . The shear modulus 
/X ijstrohmaver et al. Ill99ll ) is 



0.1194 



i + o.595(ro/r)2 



(2) 



^ The recent calculation by iLevin I ll2007h is particularly inter- 
esting, as it examines the time-dependence of the interaction be- 
tween the crust and the core in response to an initial perturba- 
tion. The behaviour is complicated, with time-dependence that is 
not pure oscillatory, but results in periodic amplification at the 
natural crust frequencies. 



Here Z is the atomic number of the ions, rii the density 
of ions and a = {S/ATvrii)^^''^ the average inter-ion spacing. 
The parameter F — {ZeY / aksT , where T is the temper- 
ature and F = Fo = 173 marks the point at which the 
solid lattice melts to form an ocean (jFarouki fc Hamaguchi I 
[1993), and we use this to determine the upper boundary of 
the crust. Crust temperature in the tail of a giant flare is 
not well constrained: o bservations set a lower limit of lO'' 
K (|Tiengo et al. I [2OO5I ') . but theory suggests that it could 
be as high as ~ 10^ K (jLvubarskv. Eichler fc Thompsoru 
[2002' ). so we examine the range T = Iff - 10^ K. 

We consider stellar models with M = 1.2 — 2.5Mq 
and R = 8 — 15 km, subject to causality constraints, 
to cover t he full range of possi ble strange star parame- 
ter space l|Page fc Reddv I I2OO6I ). Both thin nuclear and 
nugg e t crust models are studied (lAlcock. Farhi fc Olinto I 



Il98d : IJaikumar. Reddv fc Steiner 1120061 ) We model mag 
netic fields in the range B = 10^^ - lO^'' G. The strong 
field a ffects only the very ou termost low density layers of the 
crust (|Harding fc Laill2006h : we estimate the corrections to 
computed mode frequencies to be less than 1%. Energetic 
arguments show that global crust structure, including the 
suspended crust model, is unaffected. 

For the thin nuclear crust there is some uncertainty in 
the composition of neutron-rich nuclei, which will be con- 
strained by future rare isotope experi ments. We therefore 
survey a range of equations of state (iHaensel fc Pichon I 



ll994l : iRuster, Hempel fc Schaffner-Bielich 1 12006''). 

The nugget crust is composed of a lattice of strange 
quark nuggets embedded in a background degenerate elec- 
tron gas. The electrons contribute to the pressure, and the 
nuggets to the energy density. The density within the crust is 
given by p = a;eo, where eo = nquarkPg is the energy density 
of quark matter inside nuggets and depends on the density 
JT-quark ~ 1 fm^^ aud chcmlcal potential fig ~ 300 MeV at 
which stable quark matter vanishes, x is the volume fraction 
occupied by nuggets and is given by 



(3) 



where fie is the electron chemical potential, uq is the electric 
charge dens ity of the quark nugget and yq is it s charge sus- 
ceptibility (|jaikumar. Reddv fc Steiner |2006|) . The quark 
matter parameters are poorly known and can only be de- 
termined within the context of specific models. In the Bag 
model, ng = m'lfiq/2-K^ and XQ ~ A'g/^^- Further, requir- 
ing SQM to be absolutely stable and simultaneously requir- 
ing normal nuclei to be metastable restricts fiq to a narrow 
range centered around fiq ^ 300 MeV. Thus within the Bag 
model, the remaining uncertainty is parametrized through 
the effective strange quark mass rris which we expect to be 
in the range 150 - 250 MeV. 

The spherical nugget phase occupies most of the crust; 
we will neglect the small region at the base of the crust 
containing the pasta phase (see below). In computing shear 
modulus (Eqn. [2| the quantity 

Z=^fRl{nQ-XQi^.) (4) 

is the charge of the nugget where Rn = y Xd is the typical 
nugget size, = l/(47raxQ)^''^ being the Debye screen- 
ing length and a the fine structure constant. The factor 
/ is a correction due to screening inside nuggets, / « 
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Figure 1. Shear and Alfven velocities (vs, va) in the crust for 
stars with M = IAMq, R = 12 km, T = 10* K and B = 10^* G. 
Top: Thin nuclear crust (crust thickness AR/R = 5.4%). Bottom: 
Crust with nuggets, for = 150 MeV {AR/R = 0.7%) and 
rris = 250 MeV (AR/R = 6.1%). 



3(2/ - tanhy)/y^, where y « 1.6 (|Alford et al. 1 120061 ') gives 
/ ~ 0.5 for the larger droplets. The quantity a = R„/x^^^ 
is the average inter-nugget distance, and rii = 3/(47ra^) the 
density of nuggets. 

The variation of with depth is given by integrating 
the equations of h ydrostatic equilibrium. Using the l imits 
a:: = and x — 1. Ijaikumar. Reddv fc Steiner I (|2006t ) fol- 
lowed this procedure to estimate A_R. In reality the region 
with X > 0.5 that contains the pasta phase occupies a tiny 
region with thickness < Im. Further, with increasing x the 
free energy gain of the heterogeneous state becomes negligi- 
ble and even a small surface tension can disfavor the large 
x region. For these reasons we set x ~ 0.5 at the base of 
the crust. The value of x at the top of the crust is set by 
the melting condition: for the parameters examined, it lies 
in the range x - 10"^^ - 10"*. 

Figure[T]shows the variation of Vs and va with depth f or 
example crust models (compare to Fig. 1 of IPiro I l|2005h ). 
The shear speed in the nugget crust is smaller than in the 
thin-nuclear crust. This can be understood by noting that 
at constant pressure Vs ~ \/ Z^/^/A where A denotes the 
baryon number. Further, and unlike in the nuclear case, both 
Z and Z/A of the nuggets decrease rapidly with depth. 



2 RESULTS 

Our results are good agreement with a n alytic est i mates 
jMcDermott. van Horn &: Hansen! 1 19881: lOuncan I 1 19981 : 



IPiro I I2OO5I : ISamuelsson fc Anderssonll2007l ). Overall, fre- 
quencies scale with Vs- The n = mode frequencies are 
almost independent of AR, and varying B and T changes 
frequency by no more than a few Hz. The frequencies scale 
as [{l + 2){l — l)]^^'^, making it difficult to fit a mode sequence 
with an 18 or 26 Hz fundamental. Following previous stud- 
ies we assume that these represent global Alfven modes, and 
search for a fundamental at « 30 Hz. The radial overtones 
vary little with I, but depend strongly on AR, frequencies 
tending to increase as the crust thins. Crust thickness is set 
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Figure 2. Thin nuclear crust mode frequencies (for a star with 
M = IAMq, R = 12 km). Left: effect of varying magnetic field, 
for fixed T = 10* K. Right: effect of varying temperature, for fixed 
B = 10^2 Q Crust thickness AR/R falls from 5.4% at T = 10* 
K to 2.1% at T = 10^ K. 



in part by compactness. Increasing temperature also tends 
to push overtone frequencies up as the outer layers of the 
crust melt. The overtones also depend strongly on B, since 
Va exceeds Vs at magnetar field strengths. 

For the thin nuclear crust models the residual uncer- 
tainty in the equation of state leads to variations of at most 
a few percent, and i n wha t follows we quote results based on 
iHaensel fc Pichon~l (|l994l '). For the parameter space studied 
the frequency of the n = 0, / = 2 mode lies in the range 
26 - 54 Hz. To obtain a fundamental ^ 30 Hz requires 
high mass (M ^ 2.4M0 for aU radii, or M ^ 2.2Mq and 
R ^ 1 4 km), pushing the li mits of strange star parameter 
space (jPage fc Reddv II2OO6I ). 

The overtone frequencies are all high. Even at B = 10^^ 
G, the lowest frequency is « 725 Hz, for the model with the 
thickest crust (M = I.2M0, 7? = 15 km, AR/R = 8%). For 
B > 10^* G, the lower limit on overtone frequency is higher 
still, at « 1100 Hz. There is no model that would allow an 
overtone at 625 Hz. Figure [2] illustrates the effects of varying 
B and T on the frequencies of both the n = and n = 1 
modes. 

For the nugget crust models there is additional uncer- 
tainty in the value of tUs. Increasing rria increases shear 
speed, and for the n — modes frequency scales directly 
with TUs. However, in general Vs is lower than for the nuclear 
crust models. For the parameter space studied the frequency 
of the n = 0,1 — 2 fundamental lies in the range ~ 1 — 21 
Hz. There is no model that would permit a fundamental in 
the range 28-30 Hz. 

The range of possible frequencies for the n = 1 over- 
tones is extremely large. Crust thickness increases dramat- 
ically with nis, and this effect can in fact compensate for 
the change in Vs- For the highest value of nis, AR/R is 
comparable to or thicker than the thin nuclear crust; for 
the lowest values, however, it can be an order of magni- 
tude thinner. Temperature dependence is complex. As T in- 
creases from lO'^ K frequency drops as Vs drops, the effect 
being more pronounced than for the nuclear crust. However 
at high enough T crust thinning accelerates dramatically. 
This can offset the effects on Vs and cause frequency to rise 
again. At B = 10^^ G, frequencies for the n = 1 overtone 
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Figure 3. Nugget crust mode frequencies (for a star with M = 
IAMq, R = 12 km) for strange quark mass ma = 150 MeV (solid) 
and rus = 250 MeV (dotted). Left; effect of varying magnetic 
field, for fixed T = 10* K. Right: effect of varying temperature, 
for fixed B = 10^^ G. The minimum in frequency for the overtone 
occurs at higher temperatures for rris = 250 MeV. 

can be as low as ~ 260 Hz. Frequencies do however rise 
as B increases, and by B = 10^* G only a few very ex- 
treme models permit an overtone frequency as low as 625 
Hz {M = 1.2M0, i? = 15 km, tUs ~ 250 MeV, T ^ 9 x 10* 
K). Figure [3] illustrates the effects of varying B, T and tUs 
on the n = and n = 1 modes. 



3 DISCUSSION 

For neutron star models torsional shear modes, or their 
global magneto-elastic equivalents, are a good fit for the 
observations. For strange stars, the situation is much more 
difficult. A thin nuclear crust model can give a fundamental 
in the right range if stellar mass is large, but the overtone 
frequencies are far too high. The nugget crust model per- 
mits a wider range of frequencies due to the uncertainty in 
the strange quark mass mg. The lowest order n = modes 
could explain some of the QPOs in the range 18-150 Hz, but 
because the fundamental frequencies are so low it is difficult 
to fit a mode sequence, given t he expected scaling with I 
ISamuelsson fc Anderssonl l2007l ) . For the radial overtones 
there are regions of parameter space at high temperature 
with modes at the right frequency, but only for magnetic 
fields at least an order of magnitude lower than those in- 
ferred for magnetars. At magnetar field strengths only mod- 
els at the very limits of parameter space permit an overtone 
frequency as low as 625 Hz. The other constraint on the 
nugget crust is the sensitivity of the overtones to temper- 
ature fiuctuations (greater for this model than for the thin 
nuclear crust or neutron star crusts). The observations in- 
dicate that the 625 Hz QPO has high coherence and lasts 
for several hundred seconds, a period duri ng which tcmper- 
ature could vary substantia lly (|Watts fc S trohmavcr 200^; 
IStrohmaver fc Watts II2OO6I ). 

We conclude that the frequencies of toroidal 
shear modes in strange star crusts have seri- 
ous difficulty explaining the QPO frequencies ob- 
served during magnetar hyperflares. I f the res ults of 
iGlampedakis. Samuelsson fc Andersson I l|2006l ) and iLevin I 
( 20071 ) hold true, this conclusion will not change greatly 



when coupling to core is included in the calculation. 
However, these results must now be verified using more 
sophisticated models that include a realistic stellar geome- 
try, field configuration, and general relativistic corrections. 
There may also be other types of modes in the right 
frequency range l|Chugunov 1 1200^ ) . although these should 
be harder to excite and detect. We note that there are 
alternative non-seismic models for the QPOs, but that 
these models have serious d ifficulties (for a discussion see 
IWatts fc Strohmaver I (|2007l )'). 

The clear distinction between the theoretical predic- 
tions for neutron star and strange star crust models is ex- 
tremely promising. It offers a robust, largely model inde- 
pendent means of distinguishing strange stars from neutron 
stars, something that has in the past been lacking. This type 
of study should lead to rapid progress in constraining the 
equation of state of compact stars. Ruling out the strange 
star hypothesis would directly impact the phase diagram of 
QCD at finite chemical potential. It would offer a strong 
constraint on models of dense quark matter, indicating that 
the deconfinement transition is not significantly lowered by 
the dynamics of the strange quarks. One important impli- 
cation is that multiply strange hadronic states such as the 
H-dibaryon in the terrestrial context are less likely. In ad- 
dition, constraints on the models at finite chemical poten- 
tial should be relevant for finite temperature extensions of 
these models which have been employed to describe heavy- 
ion experiments such as the Relativistic Heavy-Ion Collider 
(RHIC). Perhaps most importantly, ruling out strange stars 
directly from observations is the only way to ascertain that 
terrestrial matter is stable. 
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